Lebed-PRL-2008 



Paramagnetic Intrinsic Meissner Effect in Layered 

Superconductors 

A.G. Lebed* 
Department of Physics, University of Arizona, 
1118 E. 4-th Street, Tucson, AZ 85721, USA 

Abstract 

Free energy of a layered superconductor with £j_ < d is calculated in a parallel magnetic field 
by means of the Gor'kov equations, where £j_ is a coherence length perpendicular to the layers 
and d is an inter-layer distance. The free energy is shown to differ from that in the textbook 
Lawrence-Doniach model at high fields, where the Meissner currents are found to create an 
unexpected positive magnetic moment due to shrinking of the Cooper pairs "sizes" by a magnetic 
field. This paramagnetic intrinsic Meissner effect in a bulk is suggested to detect by measuring 
in-plane torque, the upper critical field, and magnetization in layered organic and high-T c 
superconductors as well as in superconducting superlattices. 
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The Meissner diamagnetic effect is known to be the most important property of a super- 
conducting phase and is responsible for destruction of superconductivity both in type-I and 
type-II superconductors [1]. Meanwhile, as shown by us in Refs. [2-4] and independently 
by Tesanovic, Rasolt, and Xing in Ref. [5], quantum effects of an electron motion in a 
magnetic field result in the appearance of a qualitatively different phenomenon - supercon- 
ductivity surviving in high magnetic fields in layered [2-4] and isotropic three-dimensional 
[5] type-II superconductors. In particular, it was shown [2-4,6] that in a layered conduc- 
tor in a parallel magnetic field, where the Landau levels quantization is impossible, some 
other quantum effects - the Bragg reflections - play an important role. These quantum 
effects result in a "two-dimensionalization" (i.e., 3D — > 2D dimensional crossover) of an 
open electron spectrum in an arbitrary weak parallel magnetic field. This is known [6,7] to 
cause the field-induced spin- density- wave (FISDW) and field-induced charge-density-wave 
(FICDW) instabilities in layered quasi-one-dimensional (Q1D) conductors. More compli- 
cated 3D — > ID — > 2D dimensional crossovers are shown [8,9] to be responsible for the 
experimentally observed non-trivial angular magnetic oscillations in a metallic phase of dif- 
ferent layered organic conductors, including Interference Commensurate [8] and Lebed Magic 
Angle [9] oscillations. 

As shown in Refs. [2-4,7], the similar 3D — > 2D dimensional crossovers have to be re- 
sponsible for a stabilization of a superconducting phase in layered Q1D [2,4] and quasi-two- 
dimensional (Q2D) [3] conductors since 2D superconductivity is not destroyed in a parallel 
magnetic field. More precisely, it is shown [2-4] that: (i) the quantum effects make the 
upper critical field to be divergent, H^{T) — > oo as T — > 0, and (ii) there is some critical 
filed, H*, above which superconducting temperature grows in an increasing magnetic field. 
Such superconducting phase with dT c /dH > is called the Reentrant Superconductivity 
(RS) [2-5]. The original predictions [2-4] have been theoretically confirmed by a number of 
studies [10-14], including a study [14], which takes into account a possibility of a pure 2D 
phase transition. Despite of a great success of 3D — > ID — > 2D and 3D — > 2D dimensional 
crossovers concepts in the explanations of magnetic properties in a metallic [7-9] , the FISDW 
[6,7], and the FICDW [7] phases of organic conductors, so far there has been no evidence 
that superconducting temperature can grow in high magnetic fields due to the the quantum 
3D — > 2D dimensional crossovers [2-5,10-14]. A possibility of the RS phase to exist was 
experimentally studied in Q1D layered organic superconductors (TMTSF) 2 X (X = PF 6 and 
X = C10 4 ) by Naughton's and Chakin's groups [15-17]. Their experiments gave hints on 
a possibility for superconductivity to exceed significantly the quasi-classical upper critical 
field -ffj 2 (0) - the effect predicted in Refs. [2-5,10-14] - but they were not able to confirm 
the appearance of the RS phase with dT c /dH > 0. Analogous experiments, performed on a 
Q2D superconductor Sr 2 Ru04 [18], did not detect any stabilization of superconductivity at 
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The main obvious difficulty in the above mentioned efforts to discover the RS phase is 
the Pauli spin-splitting destructive mechanism against superconductivity and the related 
Clogston paramagnetic limiting field, H p [1]. It is absent only for some triplet supercon- 
ducting phases which are believed to exist in (TMTSF) 2 X [2,15,4] and Sr 2 Ru0 4 [19] super- 
conductors. On the other hand, recently there have appeared the NMR measurements [20] 
in favor of a singlet nature of superconductivity in (TMTSF) 2 C10 4 material as well as some 
doubts [21] in a triplet nature of superconductivity in Sr 2 Ru04 one. 

The goal of our Letter is a three-fold one. Firstly, we show that, although in layered 
paramagnetically limited (singlet) superconductors the RS phase may not be characterized 
by dT c /dH > feature [2-4, 10-14], nevertheless the RS phase reveals itself as another 
unique phenomenon - paramagnetic intrinsic Meissner effect (PIME). Secondly, we extend 
microscopical theory [3] to describe the most important from experimental point of view 
d- and s-wave Q2D superconductors with £j_ < d, where £j_ is a coherence length perpen- 
dicular to the conducting layers and d is an inter-layer distance. And thirdly, we suggest 
simple experimental methods to detect the PIME phenomenon in Q2D organic and high- 
T c superconductors by using in-plane torque, the upper critical field, and magnetization 
measurements. In particular, we demonstrate that in-plane anisotropy due to anisotropic 
Ginzburg-Landau coherence lengths, which disappears in an intermediate region of magnetic 
fields (where the Lawrence-Doniach model is applicable), appears again in high magnetic 
fields as a consequence of the PIME phenomenon (see Figs. 1,2). We suggest to measure 
in-plane anisotropy of the upper critical field and magnetization as well as in-plane torque 
in high magnetic fields to discover the PIME and RS phenomena. For these purposes, we 
derive a free energy of a Q2D superconductor with £j_ < d in a parallel magnetic field from 
the Gor'kov formulation [22-24] of the microscopic superconductivity theory. Our results 
coincide with that of the Lawrence-Doniach model [25,26] only at low enough magnetic 
fields, H <C H*, where the Meissner effect is diamagnetic. We show that, at high magnetic 
fields, H ~ H* < H p , the field starts to shrink the Cooper pair "sizes" in perpendicular to 
conducting layers direction due to 3-D — > 2D crossovers in a parallel magnetic field. The 
above mentioned 3-D — > 2D crossovers of the Cooper pairs are not taken into account in 
the Lawrence-Doniach model and, as shown below, are responsible for the unique PIME 
phenomenon. 

Let us consider a layered superconductor with a Q2D electron spectrum, 



e(p) = H (Px, V y ) + 2t_L cos(p z d) , t± < e F , 



(1) 



in a parallel magnetic field, 



H = (0,#,0), A = (0, 0, —Hx) , 



(2) 
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where €\\(p x ,Py) ~ £f is in-plane electron energy, t± is an overlapping integral of electron 
wave functions in a perpendicular to the conducting planes direction, e F is the Fermi energy. 
Electron spectrum (1) can be linearized near 2D Fermi surface (FS), e\\(p x ,p y ) = e F) in the 
following way, 

e(p) -e F = v x (p y )\p x - p x (p y )] + 2t± cos(p z d) , (3) 

where v x {p y ) = de\\(p x ,p y )/dp y is a velocity component and p x {p y ) is the Fermi momentum 
along a;-axis. 

In the gauge (2), electron Hamiltonian in a magnetic field can be obtained from Eq.(3) by 
means of the Peierls substitution method, p x — > —i(d/dx), p z — > p z + (e/c)Hx [6]. Therefore, 
electron Green functions in a magnetic field satisfy the following differential equation, 



iu n - V x (p y 



+2t ± cos[p z d+——)+2^ B Hsj 

xG iu!n (x,xi,p y ,p z ] s) = S(x - xi) , 



(4) 



where u n is the Matsubara frequency [22], /i B is the Bohr magneton, s = ±| is an electron 
spin projection along quantization y-axis; H = 1. It is important that Eq.(4) can be solved 
analytically. As a result, we obtain, 

co n (x - Xi)~ 



Gi^x.x^py.p^s) = -% — N exp 



x exp 



2ifi B sH(x — xi) 

V X (Py) 



exp 



V x (Py) I V X (Py) 

iHPv) r • ( i eHdx\ . ( , 

~m{p z d-\ \ — sm\p z d + 



exp[ip x (p y )(x - x x )\ 
eHdxi " 



(5) 



where X(p y ) = 4t±c/ev x (p y )Hd. 

Linearized gap equation determining superconducting transition temperature, T C (H), can 
be derived using Gor'kov equations for non-uniform superconductivity [3,23,24]. As a result, 
we obtain, 



A(x) = V 



dl 



dxi 



2ttT 



v ± (l) J\ x - xl \>\ Vx (i)\/n ^(/) sinh 



27rT|a-xi| 
V x (l) 



COS 



2fj, B H(x - x ± ) 



v x (l) 



x J < 2A(Z) sin 



eHd(x — xi) 
2~c 



sin 



eHd(x + xi) 



2c 



A(x x ) , 



(6) 



where integration in Eq.(6) is made along 2D contour, €\\(p x ,p y ) = e F , v±(l) is a veloc- 
ity component perpendicular to the 2D FS, V is an effective electron-electron interactions 
constant, Q is a cut-off energy. [Note that, although Eq.(6) is derived for singlet s-wave 
superconductors, it is also valid for <i-wave superconductors [27] if we redefine properly 
anisotropic coherence lengths and the effective interactions constant V '.] 
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We point out that Eq.(6) is the most general one among the existing equations to de- 
termine the parallel upper critical field in a layered superconductor. In particular, it takes 
into account the Bragg reflections and related 3-D — > 2D dimensional crossovers of elec- 
trons, which move in the extended Brillouin zone in a parallel magnetic field. As shown 
in Ref. [7], the above mentioned quantum effects result in a momentum quantization law 
for an electron momentum component along rr-axis. This a reason why the kernel of the 
integral Eq.(6) is periodic [2,4] with respect to variables x and x\. In the case, where the 
destructive Pauli spin-splitting effects against superconductivity are absent [i.e., at hb — 
in Eq.(6)], Eq.(6) possesses a periodic solution for A(x) at any magnetic field. In this case, 
superconductivity is stable in an arbitrary strong magnetic field and exists at high fields 
in a form of the RS phase with dT c /dH > 0. In the case of a singlet superconductivity, 
which is considered in the Letter, the Pauli spin-splitting effects may eliminate the super- 
conductivity with dT c /dH > 0. Nevertheles, in the latter case, the RS phase reveals itself 
as unusual anisotropy of the upper critical field and magnetization in high magnetic fields, 
H>H*~ (tJ^yPfo/Ccd < H p (see Figs. 1,2). 

As the most general equation, Eq.(6) contains Ginzburg-Landau and Lawrence- Doniach 
descriptions as its limiting cases at low enough magnetic fields, H <C H*. For the so-called 
Josephson coupled layered superconductors with £j_ < d [25,26], Eq.(6) may be simplified 
and rewritten in the following differential form, 



T r -T 



T, 



- 2.1 



TrT r 



d 



/2xuj r 



with 



and 



where 



Mm = a 4 



i 2 



V x (l) 



dz 



o sinh(z) 



sm 



u c v x (l) 



A(x) = 



4ttT c v f 



B{H) 



8*1 



vf 
v x (l) 



dz 



o sinh(,2) 
dl 



sm 



u c v x (l) 



4ttT c v f 
dl 



cos 



AixT c v F 



(7) 
(8) 
(9) 

v ± (l)VJ/ J v ± (l) ■ (10) 
[Here, u c = eHv F d/c is a characteristic frequency of an electron motion along open FS 
(1) [3], £ x = y / 7C(3)(f^(/)) 1 ^ 2 /47rT c is in-plane Ginzburg-Landau coherence length, fi B H ~ 
w c (H) < ttT c .] 

Note that Eqs.(7)-(10) extend the Lawrence- Doniach model [25,26] to the case of strong 
magnetic fields and can be called extended Lawrence-Doniach equations. In contrast to 
the traditional Lawrence-Doniach equations, the coefficients A(H) and B{H) in Eqs.(7)- 
(10) depend on a magnetic field, which means that a probability for the Cooper pair to 
jump from one conducting layer to another depends on the field. This important feature 
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of Eqs.(7)-(10) is a consequence of shrinking of the Cooper pairs "sizes" due to 3D — > 2D 
dimensional crossover in a parallel magnetic field [2,3,7]. 

Below, we are interested in descriptions of the RS and PIME phenomena, therefore, we 
consider Eqs.(7)-(10) at high magnetic fields. It is possible to show that at H > H* the 
solution of Eq.(7) can be represented as A(x) = A = const, which corresponds to the RS 
phase [2,3]. In this case, the corresponding second order term of a free energy with respect 
to the order parameter A can be written in the following simple form, 

T C (H)-T] 



F 2 (T,H) = -N(e F ) 

where 



T 

J- c 



(11) 



Tcm . Tc _, 1 (^_, 1 _^ + 0,5^(£f) 2 e , , 12 , 

where N{ef) is a density of states per one electron spin projection at e = e F - 

Note that the the first term in Eq.(12) describes destruction of a singlet superconductivity 
by the Pauli spin-splitting effects, whereas the last term in Eq.(12) is responsible for the 
restoration of superconductivity at high magnetic fields and for the appearance of the RS 
phase and PIME phenomenon. If we take into account that the fourth order term of a free 
energy with respect to the order parameter A can be calculated at H > H* in a standard 
manner [23], F 4 = 7C(3)iV(e_F) A 4 /167r 2 T c 2 , then we can minimize the total free energy and 
find that 

4vr 2 

F(T, H) = -j^jN(e F ) \T C {H) - T] 2 . (13) 

Magnetization can be found by a differentiation of the free energy (13) with respect to a 
magnetic field, 



8F(T,H) 8 ,[T C -T\l ln S , / efxjj, ^ 



M(T,H) = —^ = m N ((F 



-4.2 n% + 1.9 



H , (14) 



Eqs.(12)-(14) are the main results of the Letter. Note that, in Eq.(14), the first term 
corresponds to a destruction of superconductivity due to the Pauli spin-splitting effects, 
whereas the second term represents unusual paramagnetic orbital contribution to a magnetic 
moment (i.e., the PIME phenomenon). It is important that £ x in Eqs.(12)-(14) is anisotropic 
and depends on a direction of a magnetic field, since it is in-plane component of a coherence 
length perpendicular to the field. Therefore, the RS and PIME effects in Eqs.(12)-(14) can 
be detected by measuring a torque provided that spin-splitting effects are isotropic (see 
Eqs.1,2). 

In conclusion, we discuss possible experiments to discover the PIME and RS phenom- 
ena. The most direct method is to create such layered superconducting super-lattice, where 
uj c (H) 3> [jlbH [28]. The latter condition means that the orbital effects are more important 
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than the Pauli spin-splitting ones. Therefore, in this case, the increase of transition temper- 
ature (12) and the paramagnetic Meissner effect (14) can be directly observed. Nevertheless, 
in most real physical compounds with £j_ < d, uj c (H) ~ fi B H and, thus, the PIME (14) and 
RS (12) phenomena can be observed only indirectly - by measurements of anisotropies of 
the in-plane upper critical field (12) and magnetization (14) as well as by measurements 
of in-plane torque. In our opinion, the most perspective superconductors for indirect ob- 
servations of the PIME phenomenon in steady magnetic fields are organic compounds a- 
(ET) 2 NH 4 Hg(SCN) 4 , /s>(ET) 2 Cu(NCS) 2 , k-(ET) 2 Cu[N(CN) 2 ]X, a-(ET) 2 KHg(SCN) 4 , and 
A-(BETS) 2 FeCl 4 [29]. The above mentioned studies of the in-plane anisotropies can be per- 
formed also in high-temperature superconductor YiBa 2 C3C>7 but it will require ultra-high 
pulsed magnetic fields. 

The author is thankful to N.N. Bagmet, P.M. Chaikin, and M.V. Kartsovnik for the 
useful discussions. 
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FIG. 1: Superconducting transition temperature in a parallel magnetic field for a paramagnetically 
limited Q2D superconductor is sketched. GL - area of applicability of the Ginzburg-Landau theory 
[1], LD - area of applicability of the Lawrence-Doniach model [25,26], PIME - area, where both the 
GL and LD descriptions are broken. In the latter case, which corresponds to shrinking of the Cooper 
pairs "sizes" by a magnetic field, our Eqs. (6)-(14) are still valid and the Reentrant Superconducting 
(RS) phase appears. The RS phase may reveal itself as an increase of the transition temperature 
in a magnetic field, if the orbital effects of an electron motion are stronger than the Pauli spin- 
splitting effects (dashed line). The RS phase always reveals itself as a paramagnetic intrinsic 
Meissner effect (PIME), which results in unexpected in-plane anisotropy of the upper critical field 
and magnetization even in the case, where the Pauli spin-splitting effects are strong and, thus, the 
area with dT c /dH > is absent (solid line). We suggest to measure in-plane torque, the upper 
critical field, and magnetization to discover the RS phase. 
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FIG. 2: Solid line: in-plane magnetization, 47rM, is sketched as a function of a magnetic field in 
the absence of the Pauli spin-splitting effects. At high magnetic fields, H > H*, the Reentrant 
Superconducting (RS) phase reveals itself as a paramagnetic intrinsic Meissner effect (PIME). 
Dashed line: an absolute value of in-plane torque, |r|, is sketched. It is important that the torque 
is independent on the Pauli spin-splitting effects since they are isotropic. Therefore, even in the 
case, where the destructive Pauli spin-splitting effects eliminate a positive sign of the Meissner 
effect in high magnetic field, the PIME phenomenon and the RS phase can still be detected by the 
in-plane torque measurements. 



10 



